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Abstract 

Novel Positron Emission Tomography system, based on plastic scin¬ 
tillators, is developed by the J-PET collaboration. In order to optimize 
geometrical configuration of built device, advanced computer simulations 
are performed. Detailed study is presented of background given by acci¬ 
dental coincidences and multiple scattering of gamma quanta. 


1 Introduction 

The GEANT4 Application for Tomographic Emission (GATE [T]) represents 
one of the most advanced specialized software packages for simulations of PET 
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scanners. Despite the complexity of the simulated system, GATE is easily con¬ 
figurable and facilitates convenient use of the powerful GEANT4 simulation 
toolkit. 

Thanks to the fact, that the software was widely verified, it may be used 
for simulations of such a prototype devices as Strip-PET scanner [2]-Il]: build 
by the J-PET collaboration. The scanner is based on the plastic scintillators 
representing innovative approach in the field of PET tomography. Another im¬ 
portant feature of the scanner is large axial field-of-view (AEOV). PET scanners 
with large AEOV are also developed by other collaborations M-m- 

2 Setting parameteres of the simulations in the 
GATE software 

Properties of the scintillating material and the detecting surface, were set using 
three GATE-specific files: GateMaterials.dh^ Materials.xml and Surfaces.xml. 
Some of them, could be fixed using data from documentation prepared by the 
producers of the equipment. 

Eor example the properties of the scintillating material EJ230 m, that is 
used by the collaboration in real-life experiments are: 

• scintillation yield - 9,700 1/MeV 

• refraction index - 1.58 

• density 1.023 g/cm^ 

• emission spectrum - Eig. maximum of emission at 391 nm 
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Eigure 1: Emission spectrum of the EJ230 material m and quantum efficiency 
of the R4998 photomultiplier [17] 

The only property of the detecting surface (which immitates the photomul¬ 
tiplier Hamamatsu R4998 US]), that has to be set by the user is the dependence 
of quantum efficiency on the wavelength of optical photons [m (Fig. 0- 

Some important properties, however, are not given by the producer. One of 
them is the absorption length dependence on the light wavelength. Therefore, we 
adopted and tested this dependence from another similar material as described 
below. 
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2.1 Simulations of the single strip 

The dependence of absorption length on light wavelength for plastic scintillator 
may be found in Ref. [18]. Similar borrowing has been applied by the authors 
simulating the NEMO detector [19]. The dependence taken from publication was 
read out from the picture, smoothed using line interpolation and implemented in 
GATE software. After that some simulations of the single strip were performed 
and their results were compared with the experiment. 

In the experiment, the collimated source of gamma quanta (Na-22) was 
moved along the scintillator EJ230 (5 mm x 19 mm x 30 cm) with step 3 mm and 
the beam was directed perpendicularly into the scintillator. Eor each position of 
the beam, 1-dimensional histogram of the number of photoelectrons was created. 
The histogram was put into the single column of two-dimensional histogram 
presented in the background of the Eig. It- In this figure, one can see the 
dependence between the number of photons detected by the photomultiplier and 
the position of the beam of gamma quanta. Experimental data are available for 
positions between -14.7 cm and 14.7 cm and the width of bins is 3 mm. In this 
figure results for R4998 photomultiplier attached to the scintillator at the end 
(position 15 cm) are shown. 
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Eigure 2: The left part of the figure presents the dependence of the absorp¬ 
tion length on the light wavelength. Upper line shows result obtained for pure 
polystyrene (PST) [18] and the lower line is scaled by factor of 0.55. The right 
part of the figure shows the spectra of number of photoelectrons as a position 
of the beam of gamma quanta. Dashed line indicates maximum number of pho¬ 
toelectrons produced by 511 keV gamma quanta as a function of position of 
irradiation assuming absorption as measured for PST (dashed line) and PST 
absorption scaled by factor of 0.55 (solid line). 


Fig. |3] shows comparison of simulated and experimental distributions of 
number of photoelectrons for three exemplary positions. A good agreement was 
obtained when scaling the absorption length of pure polystyrene [18] by factor 
of 0.55 (Eig. [^. The scaling factor accounts effectively for the absorption due 
to the primary and secondary admixture in the scintillator material, imperfec¬ 
tions of surfaces and reflectivity of the foil. Dashed and solid line in right part 
of Eig. presents results of simulations performed for energy loss of 341 keV 
corresponding to the maximum energy of the electron scattered by the 511 keV 
gamma quantum via Compton effect. Dashed line was obtained assuming the 
absorption length as determined for the pure PST, whereas solid line shows 
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Figure 3: Comparison of the simulated and experimental histograms of energy 
deposited by 511 keV gamma quanta (in number of photoelectrons) for the beam 
positions -12 cm (left), 0 cm (middle), 12 cm (right); experimental spectrum is 
suppressed at low values due to the triggering conditions [2] . 


result after scaling the absorption by a factor of 0.55. The scaling factor was 
optimised to the experimental results. 

3 Simulations of the single layer J-PET scanner 

A diagnostic chamber of the J-PET detector will form a cylinder which will be 
constructed from the plastic scintillator strips [20]-[22]. In this article we present 
simulations for the detector with the inner radius of R=427.8 mm (radius similar 
to commercially available PET systems m, m)- We assume that the detector 
possesses one layer build out of 384 EJ230 scintillator strips with dimensions of 
7 mm X 19 mm x L (L = 20 cm, 50 cm, 100 cm or 200 cm). Geometry of the 
simulated scanner is visualised in Eig. [^ 



Eigure 4: Visualisation of the geometry of the single-layer J-PET scanner with 
radius of the cylinder R and the length of the scintillators L. 


3.1 Scattered coincidences 

In order to estimate secondary scattering of gamma quanta in the detector 
material, we have simulated annihilations homogeneously in the 2 m long line 
placed along the central axis of the scanner. In the following, we consider 
few most probable responses of the detector system (see Eig. [^. In the most 
probable case both gamma quanta will escape detection and no signal will be 
observed {Ngtrips = 0)- The second frequent category corresponds to events 
when only one strip was hit {Nstrips = !)• Eurther on for the multiplicity 
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of strips Nstrips >= 2 we can distinguish different cases for the same value of 
^strips • Therefore for the univocal description we introduce one more parameter 
/i. Various possibilities which may occur are listed below and depicted in Fig. 

El 

• ^strips ~ 3, /i = 3 

3 quanta in 3 different strips with two secondary scatterings 

• Nstrips — 2, /i — 2 

2 quanta in 2 different strips with one secondary scattering 

• Nstrips ~ — 0 

no gamma quanta registered 

• ^strips ~ ~ f 

interaction in only one strip 

• ^strips ~ 2, /i = 2 

2 interactions in 2 different strips 

• Nstrips 3, fJj 3 

3 scatterings in 3 different strips; 2 primary and 1 secondary scattering 

It is also possible that there are 4, 5 or even more scatterings, depending on 
the energy threshold applied to each hit. 



Figure 5: Pictorial definitions of the value of multiplicity /i used further in the 
following figures. 

Histograms of the multiplicity for three different energy thresholds (0 keV, 
100 keV and 200 keV) and for four different lengths of scintillators (20 cm, 50 cm, 
100 cm and 200 cm) are presented in the Fig. Results of the simulations show 
that if energy threshold is set to 200 keV, there are no events where number of 
hits is bigger than 2. Most of scattered coincidences (with multiplicity -2) is also 
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eliminated with this energy threshold. If the energy threshold is set to 100 keV, 
for lengths of scitnillators 100 cm and 200 cm, there would be even events with 
four scatterings, which may negatively influence the quality of reconstructed 
images. 




II 



/i 


Figure 6: Histogram of the multiplicity /i for different lengths L of the diagnostic 
chamber. Meanings of different values of the multiplicity /r are described in the 
text and defined in Fig. 


For Nstrips=‘^ and Nstrips=^ and length of scintillators equal to L = 50 cm, 
histograms of time differences between subsequent hits were calculated. These 
histograms are presented in Figs[^-[^ Right panel of these figures show dis¬ 
tribution of difference between ID of hit modules {AID) as a function of hit 
time difference Q The module ID increases monotonically with the grows of the 
azimuthal angle cp (see Fig. [^. Black lines in the two-dimensional histograms 
(Figs[7|-pT|) show the boundaries between events treated as useful coincidences 
and events treated as background coincidences due to the secondary scatterings. 
A positive value of fi (2 or 3) is assigned to events above the line, which are 
treated in further analysis as true conicidences. Whereas to events below the 
line a negative value of fi (-3 or -2) is assigned since these events include sec¬ 
ondary scattering of gamma quanta. This boundary was used to separate events 
with different multiplicities for preparation of histograms presented in Fig. 

In Fig. for energy thresholds 0 keV and 100 keV, in two-dimensional his¬ 
tograms there is longitudinal structure extending between points (0 ns, 0) and 
(3 ns, 192). These events correspond to difference between time of primary re¬ 
action of the gamma quantum in a given scintillator and a time of the secondary 

IDs of hit modules are IDi and ID 2 then A/D = min{\IDi — /D 2 I, 384 — \IDi — /D 2 1). 
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scattering. The larger is the angle of the primary scattered gamma quantum the 
larger will be the AID value and also a At. For example bin with coordinates 
(2.9 ns, 192) corresponds to the backscattering - primary particle is backscat- 
tered and it is registered in the strip on the opposite side of the scintillator (2.9 
ns is the time needed by the gamma quanta to travel between opposite strips 
with speed of light). 

If the energy threshold is set to 200 keV, nearly all scattered coincidences are 
eliminated. In the lower panel of this figure there are results for this threshold. 
In ideal situation, time difference for this simulation for true coincidences would 
be always 0 and we would have only one bin for 0 ns. Because of the fact that 
gamma quanta interact with matter in different depths (Depth of Interaction), 
time difference is changing from 0 to about 80 ps. This picture show, what is 
the time limit for time-of-flight determination with scintillator strips of 19 mm 
thickness. 

In Fig. for energy thresholds 0 keV and 100 keV, in two-dimensional 
histograms there is symmetrical butterfly-shape structure extending between 
points (0 ns, 0) and (3 ns, 192) and between points (3 ns, 0) and (0 ns, 192). 
Each event with three hits and deposited energy above the energy thresold, gives 
two inputs to these histograms. An additional structure (for Nstrips = 3) which 
is spanned between points (3 ns, 0) and (0 ns, 192) originates from the time 
differences between the primary interaction of one of the gamma quantum and 
a secondary interaction of the other or from the time difference between two 
secondary interactions. Pictorial definitions of these situations are presented in 
Fig. If only the first time difference is taken into account, histograms for 
3 hits (Fig. look like histograms for 2 hits (Fig. [^. 

Response of the detector to the annihilations in the 2 m long line placed along 
the detector axis was simulated also for other lengths of scintillators L = 20 cm, 
100 cm and 200 cm. Results of these simulations for two energy thresholds 
(0 keV and 200 keV) are presented in Fig. 

scintillators, the wider the longitudal structure described above 


One can see that, the longer the 
It is caused 


by the fact that, the longer the scintillators, the longer the possible distance 
between places of the primary and secondary interactions. For the scanner 
with 20 cm scintillators, the longest possible path along the diagonal of the 
longitudinal cross-section of the scanner has length of 88 cm (^ 2.9 ns) and 
for the scanner with 200 cm scintillators, the longest possible path is equal to 
218 cm (7.3 ns). 
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Figure 7: Distributions of differences of hit times; Nstrips = 2, /r = — 2 or 
li = 2] black line in the two-dimensional histogram shows the boundary between 
events treated as originating from primary interactions only (above the line) 
and events including secondary interactions (below the line). Figure presents 
results of simulations for L=50 cm. The time differences are calculated only for 
interactions originating from the same annihilation process. 
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Figure 8: Distributions of differences of hit times; Nstrips=‘^i /i = —3 or /i = 3, 
both time differences are taken into account. Figure presents results of simu¬ 
lations for L=50 cm. The time differences are calculated only for interactions 
originating from the same annihilation process. Figure is described with details 
in the text. 
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Figure 9: Distributions of differences of hit times; A^strips=3, /i = —3or/i = 3, 
only first time difference is taken into account. Figure presents results of sim¬ 
ulations for L=50 cm. The time differences are calculated only for interactions 
originating from the same annihilation process. Figure is described with details 
in the text. 
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Figure 10: Distributions of differences of hit times (At) as a function of ID 
difference between hit modules (A/D). The length L of the detector and the 
energy thresholds are indicated above the figures. Only interactions originating 
from the same annihilation process are taken into account. 
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3.2 Accidental coincidences 

An accidental coincidence is the coincidence, in which two events occur simul¬ 
taneously in a fixed time window but in fact they are independent, they come 
from different annihilations. Because of that, number of accidental coincidences 
depends on the width of the time window, the size of the detector and in contrast 
to the secondary scattering, the accidental coincidences depend on the activity 
of the source. 


3.2.1 Accidental coincidences as a function of the source activity 


Simulations described in this section were performed for other activities of the 
source: 5 MBq, 10 MBq, 20 MBq, 30 MBq, 50 MBq, 100 MBq, 150 MBq, 
200 MBq, 250 MBq, 300 MBq, 350 MBq and 400 MBq. For each of these activ¬ 
ities 10^ annihilations were simulated. Results of simulations for the smallest 
(5 MBq) and the largest (400 MBq) activity are presented in Figs pT] and 12 
In Fig. pT| histograms contain all time differences both for hits from the same 
event and for hits from different events (there is no time window). One can see 
that the first bin is higher than expected from the general exponential depen¬ 
dence. This is because this bin contains both true and accidental coincidences. 
The structure is better visible in Fig. In the upper panel of this figure, 
histograms contain time differences between hits from the same and from differ¬ 
ent annihilations. If time differences from the same annihilations were omitted, 
there would be only accidental coincidences, as it is presented in the bottom 
panel of the figure. 
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Figure 11: Histograms of numbers of all differences of hit times for length of 
scintillators equal to 50 cm, activities of 5 MBq and 400 MBq and for two 
different energy thresholds as indicated in the legends. 
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Figure 12: Distributions of time differences between hits in the detector origi¬ 
nating from the same and from different annihilations (top panel) and from the 
different annihilations only (bottom panel). Shown are results of simulations 
for activities of 5 MBq and 400 MBq and for two different energy threshold as 
indicated in the figures. 
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3.2.2 Accidental coincidences for time windows 3 ns and 5 ns 

For simulations described in this article with the virtual linear source of annihi¬ 
lations placed along the main axis of the scanner, true coincidences are defined 
as two hits from the same annihilation having AID vs At above the black lines 
shown in Figs[^-[^ Fig. \T^ presents rate of such defined true coincidences as 
a function of annihilation source activity, time window, energy threshold, and 
detector length L. 


F,,=0keV, N^trrps='^> M=2 F,,=200keV, N^tnps='^> M=2 




Figure 13: Simulated rate of true coincidences as a function of time window, 
activity and detector length. The sequence of curves in the figure is the same as 
in the legend (from top to bottom); bottom pictures present the same data as 
the top ones but in logarithmic scale. Results for the time window of 3ns (solid 
lines) are indistinguishable from the results for time window of 5 ns (dashed 
lines). 


Accidental coincidences for time windows 3 ns and 5 ns and for four lengths 
of the scintillators are presented in Fig. One can see that if the energy 

threshold is 200 keV (right column of the figure), a rate of accidental coinci¬ 
dences is reduced by the factor of about 7 in comparison to situation when there 
is no energy threshold (left column of the figure). 

Fig. shows rate of accidental coincidences under condition that differ¬ 
ence AID is larger than 96. Which means that interactions of gamma quanta 
occurs in two different quarters of the cylinder (consisting of 384 scintillator 
strips). Such condition decrease the field of view of the detector to the cylinder 
with radius of 30 cm, however this additional condition reduces the number of 
accidental coincidences by the factor of 2. 

In Fig. rates of true and accidental coincidences are presented. The 
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Figure 14: Accidental coincidences for time windows 3 ns (solid lines) and 5 ns 
(dashed lines) for different lengths of the scintillators (as indicated in the leg¬ 
end); the sequence of curves in the figure is the same as in the legend (from 
top to bottom); bottom pictures present the same data as the top ones but in 
logarithmic scale. 


ratio is larger for longer scintillators. It is caused by the fact that for short 
scintillators there are additional accidental coincidences caused by the gamma 
quanta from outside of the tomograph. 


4 Summary 

Physical properties of the scintillating material and the photomultiplier used in 
the J-PET detector were implemented in the GATE software. The simulations 
procedures were validated by the comparison of simulated and experimental 
results for the number of photoelectron spectra. 

In previous research, studies of simplified Strip-PET scanner were presented 
[25] . Map of efficiency of 2-strip scanner was calculated and compared with the 
geometrical efficiency of such a device. In present studies, background given by 
accidental coincidences and multiple scattering of gamma quanta was investi¬ 
gated for single-layer 384-strip J-PET scanner. 

In presented simulations, the source of annihilations was assumed to be a 
2 m long line placed along the main axis of the scanner. In order to compare 
precisely obtained results with results for another devices, in the future the 
source will be simulated in accordance with NEMA-NU-2 standard m- Even 
so, it is possible to compare orders of magnitudes of calculated parameters. Eor 
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Figure 15: Accidental coincidences for time windows 3 ns (solid lines) and 5 ns 
(dashed lines) for different lengths of the scintillators (as indicated in the leg¬ 
end); minimum difference between IDs of the strips is equal to AIDth = 96; the 
sequence of curves in the figure is the same as in the legend (from top to bot¬ 
tom) ; bottom pictures present the same data as the top ones but in logarithmic 
scale. 


F,,=0keV, AID,,=96, 



E„=200key, AID„=96, N,,,,^=2 



Figure 16: Ratios between true and accidental coincidences for time window 
3 ns and minimum difference between IDs of the hit scintillator strips equal to 
96. Right figure shows results for energy threshold of Eth = 200 keV whereas 
results in left figure were obtained for Eth = 0. 


example, results obtained for 2 m long J-PET scanner for activity of 200 MBq 
are similar to these simulated for the same length RPC-PET [10]. 
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